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Abstract. So called "analogue models" use condensed matter systems (typically 
hydrodynamic) to set up an "effective metric" and to model curved-space quantum 
field theory in a physical system where all the microscopic degrees of freedom are 
well understood. Known analogue models typically lead to massless minimally 
coupled scalar fields. We present an extended "analogue space-time" programme by 
investigating a condensed-matter system — in and beyond the hydrodynamic limit 
- that is in principle capable of simulating the massive Klein-Gordon equation in 
curved spacetime. Since many elementary particles have mass, this is an essential step 
in building realistic analogue models, and an essential first step towards simulating 
quantum gravity phenomenology. Specifically, we consider the class of two-component 
BECs subject to laser-induced transitions between the components, and we show that 
this model is an example for Lorentz invariance violation due to ultraviolet physics. 
Furthermore our model suggests constraints on quantum gravity phenomenology in 
terms of the "naturalness problem" and "universality issue" . 
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1. Introduction an motivation 

The purpose of quantum gravity phenomenology (QGP) is to analyze the physical 
consequences arising from various models of quantum gravity (QG) [1, 2]. One hope for 
obtaining an experimental grasp on QG is the generic prediction arising in many (but 
not all) models that discrete space-time at the Planck scale, M P \ = 1.2 x 10 19 GeV/c 2 , 
typically induces low-energy violations of Lorentz invariance (LI). The breakdown of 
LI will manifest in a modification of the dispersion relation. We investigate Lorentz 
invariance violations (LIV) in the boost subgroup, leading to an expansion of the 
dispersion relation in momentum- dependent terms, 

E 2 = m 2 c 4 + p 2 c 2 + c 4 L M P1 p/c + m p 2 /c 2 + ^ ) , (1) 

L n>3 PI J 

where both the quantity p/(Mp\ c) and the coefficients i] n are dimensionless. $. 
The present paper focuses on non-renormalizable effective field theory (EFT) with 
ultraviolet modifications in the dispersion relation [3]. There are two aspects of this 
model theory that are interesting to look at. The so called naturalness problem is 
correlated with the appearance of r]\ and i]2 in the dispersion relation (1). The low- 
order corrections do not appear to be Planck suppressed, and would therefore always 
be dominant and in disagreement with observation. In order to avoid this problem it is 
necessary to introduce a second scaling parameter, well below the Planck scale, which 
dominates at low momentum and only affects the first and second order-terms. Such a 
scenario is not well justified within an EFT framework; in other words there is no natural 
suppression of the low-order modifications in these models. A less problematic question 
is the so called universality issue, addressing whether the LIV is particle-dependent or 
not. For a "universal" LIV the coefficients r] n are the same for all types of particles. 

In order to contribute to this debate, we have chosen a rather unconventional path: 
We investigate the energy dependent behavior of sound waves in a 2-component Bose- 
Einstein condensate. Building on the existence of so-called "analogue models" (AM) for 
minimally coupled massless fields in curved space-times [4, 5, 6, 7, 8, 9, 10], we show 
how to extend the AM to include massive particles [11], and that the dispersion relation 
(now written in terms of frequency and wavenumber) is modified by Lorentz violating 
terms at the analogue Planck scale [12]: 

uj 2 = uJq + (1 + i] 2 ) c 2 k 2 + r? 4 (JL^ 2 k' + ... . (2) 

We further calculate the dimensionless coefficients i] 2 and <q 4 for both massive and mass- 
less quasi-particles and discuss the naturalness problem and universality issue in our 
2-component AM. 

| The particular inertial frame for these dispersion relations is generally given by the Cosmological 
Microwave Background (CMB) frame. 
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2. Sound waves in 2-component BECs 

The basis for our AM is an ultra-cold dilute atomic gas of N bosons, which exist in 
two single-particle states |1) and |2). For example, we consider two different hyperfme 
states, |F = \,m F = 1) and |F = 2,m F = 2) of 87 Rb [13, 14]. They have different total 
angular momenta F and therefore slightly different energies. That permits us, from a 
theoretical point of view, to keep m\ ^ m 2 , even if they are approximately equal (to 
about one part in 10 16 ). At the assumed ultra-cold temperatures the atoms only interact 
via low-energy collisions and the 2-body atomic potential can be replaced by a contact 
potential. That leaves us with with three atom-atom coupling constants, Un, U22, and 
U12, for the interactions within and between the two hyperfme states. For our purpose 
it is essential to include an additional laser field, that drives transition between the two 
single-particle states. § The rotating frame Hamiltonian for our closed 2-component 
system is given by: || 

k = [ dr { E f-^5r^ + *1 wo^) 

J i=l,2 ^ 1 ' 

+ \ e (Uij*m*i*j + WiMijVj) } , (3) 

M=l,2 J 

with the transition energy A = hunabi containing the effective Rabi frequency between 
the two hyperfme states. Here ^(r) and ^(r) are the usual boson field annihilation 
and creation operators for a single-particle state at position r, and o x is the usual Pauli 
matrix. For temperatures at or below the critical BEC temperature, almost all atoms 
occupy the spatial modes ^i(r) and ^(r). The mean-field description for these modes, 



h 2 



V 2 + Vi-in + U ti |^| 2 + Uij 



Vi + Wj, (4) 



2mi 

are a pair of coupled Gross-Pit aevskii equations (GPE): — > (1, 2) or — > (2, 1). 

In order to use the above 2-component BEC as an AM, we have to investigate 
small perturbations (sound waves) in the condensate cloud, The excitation spectrum 
is obtained by linearizing around some background densities pio and phases 9 i0 , using: 

*i = y/pa + epn e^° +£ ^ for i = l,2. (5) 

A tedious calculation [11, 12] shows that it is convenient to introduce the following 2x2 
matrices: An effective coupling matrix, 

E = E + X, (6) 

§ A more detailed description on how to set up an external field driving the required transitions can 
be found in [15]. 

|| In general, it is possible that the collisions drive coupling to other hyperfme states. Strictly speaking 
the system is not closed, but it is legitimate to neglect this effect [16]. 

If The perturbations have to small compared to the overall size of the condensate could, so the system 
remains in equilibrium. 
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where 

3sH<( = ^^(^- ( -ir^-L.) (7, 

and 

X = [Xh = ~ k 3 ^ = = -[X],, V 2 . (8) 

2 mi 4:mip i0 

Without transitions between the two hyperfine states, where A = 0, S only contains 
the coupling constants [H]y — > Uij/h. While S is then independent of the energy of 
the perturbations, X plays a more significant role with the energy of perturbation. For 
low energetic perturbations, in the so-called hydrodynamic approximation, X can be 
neglected, X — > 0, and S — > S. 

Besides the interaction matrix, we also introduce a transition matrix, 

A S [ A ]„ = -!L*^ (-!)<«, (9) 

and a mass-density matrix, 

^ = — hdij — ■ (10) 

rrii 

The final step is to define two column vectors, 9 = [#ii,#2i] T and p = [pn,p2i] T - 
We then obtain two compact equations for the perturbation in the phases and densities, 

o= - zp-vo-ve, (ii) 

'p= - V- (D Ve + pv ) -AO, (12) 

where the background velocity v is the same in both condensates. Now combine these 
two equations into one: 

dti-- 1 6) = -d t (S^ 1 v ■ V0) - V(v 9) + V-[(D- v E' 1 v ) W] + A9. (13) 

In the next section we show how this equation is analogous to a minimally coupled scalar 
field in curved space-time. 

3. Emergent space-time in the hydrodynamic limit 

Instead of keeping the analysis general [11, 17], we now focus on the special case when S 
is independent of space and time, and on the hydrodynamic limit where X — > 0. Then 
defining 

e = ~- 1/2 o, (14) 

equation (13) simplifies to 

d 2 t 9 = -d t (l v ■ V9j —V-(v 19^+V- [(C 2 - v I vo) Vfl] + n 2 9, (15) 
where 

Cq = S 1 / 2 D S 1 / 2 ; and Q 2 = E 1/2 A S 1/2 . (16) 
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Both Cq and Q 2 are symmetric matrices. If [Cq, Q 2 ] = 0, which is equivalent to the 
matrix equation DHA = ASD, then they have common eigenvectors. Assuming (for 
the time being) simultaneous diagonalizability, decomposition onto the eigenstates of 
the system results in a pair of independent Klein-Gordon equations 

; 1 d a { y^g^i {g V u) ah d b 9 1/u \ + ul n 9yu = , (17) 
\/-gi/u L J 

where the "acoustic metrics" are given by 

r -fa-«g) i -* T ~ 



(Si/n) 



ah 



Pi/nV 

Cl/Il) 



—V I Irfxd 

and where the overall conformal factor depends on the spatial dimension d. The metric 
components depend only on the background velocity v , the background densities p i0 , 
and the speed of sound of the two eigenmodes, which is given by 

tr[C 2 ]±y/tr[C ^-4det[C 2 ] 
c i/ii — 2 ' ^ ' 

Considering the line element obtained from the acoustic metric (18), it is clear that the 

speed of sound in the AM takes the role of the speed of light. 

It is also possible to calculate the eigenfrequencies of the two phonon modes, 

J* = 0; c4 = tip 2 } . (20) 

A zero/ non-zero eigenfrequency corresponds to a zero/ non-zero mass for the phonon 
mode. They both "experience" the same space-time if 

tr[C 2 ] 2 -4det[C 2 ] =0. (21) 

The fact that we have an AM representing both massive and massless particles is 
promising for QGP if we now extend the analysis to high-energy phonon modes so 
that 1^0. For the following, we concentrate on flat Minkowski space-time, by setting 
^o = 0. 



4. QGP beyond the hydrodynamic limit 

Starting from equation (13) for a uniform condensate, we set the background velocity to 
zero, vq = 0, but keep the quantum pressure term, 1^0. The equation for the rotated 
phases 

6 = S- 1 / 2 9 (22) 
in momentum space is then [12] 

cu 2 9 = <jVs + X k 2 [D k 2 + A] VS + X A; 2 1 9 = H{k 2 ) 9. (23) 
Thus the perturbation spectrum must obey the generalized Fresnel equation: 

det{uj 2 l- H(k 2 )} = 0. (24) 
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That is, the dispersion relation for the phonon modes in a 2-component BEC is 



2 _ tr[H{k 2 )] ± y/tr[H{k 2 )] 2 -4det[H(k 2 ) 
^i/ii 2 

and a Taylor-series expansion around zero momentum gives 



co 2 - u 2 I + dUj2l/U 

^i/ii - u; i/ii| fc _ +0 + 



i d V/n 



2d(^ 



(k 2 ) 2 + O [{k 2 f] . 



(25) 



(26) 



These are two dispersion relations in the desired form of equation (2). Below we will 
explicitly compute equation (26) up to the fourth order. The fact that Uy u = Uy n (k 2 ) 
only permits even powers in k, therefore the dispersion relation is invariant under parity. 
This is by no means a surprising result, because the GPE (4) is also invariant under 
parity. 

We now define the symmetric matrices 



C 2 = Cl + AC 2 ; AC 2 = X^AX 1 / 2 ; 



Z 2 = 2X^ 2 DX 1 ' 2 = —M 



-2 



(27) 
(28) 



Note that all the relevant matrices (equations (16), (27), and (28)) have been carefully 
symmetrized, and note the important distinction between Cq and C 2 . Now define 
, 1 



(29) 



which approaches the speed of sound c 2 — * C 2 ,, in the hydrodynamic limit C 2 — > Cq, 
(see equation (19)). The second and fourth order coefficients in the dispersion relations 
(26), for a detailed calculation see [12], are: 

2tr[n 2 c 2 } - tr[n 2 } tr[c 2 } 



dk 2 
dV/n 



fc^O 



d(k 2 y 



k^O 



= C 
1 

~ 2 



1 ± 



tr[C 2 ]tr[ft 2 ] 



(30) 



2774 



± 



h 



tr[C 



tr[fi 2 ] 

2l 2 -4det[C 2 ] tr[C 2 ] 2 
tr[fi 2 ] T tr[Q 2 ] 



M, 



P i 



(31) 
(32) 



5. Lorentz violations from UV physics 



In order to obtain LIV purely due to ultraviolet physics, we demand perfect special 
relativity for X — > 0. In other words, we require all terms in the equations (30) and (32) 
which would otherwise remain in the hydrodynamic limit to be zero. The constraints 
we get are: 



Cl: tr[C 2 ] 2 -4det[C 2 ] =0; 

C2 : 2tr[fi 2 C 2 ] - tr[fi 2 ]tr[C 2 ] = . 



(33) 
(34) 
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Beyond the hydro dynamic limit, but imposing CI and C2, the equations (30) and (32) 
simplify to: 



dk 2 

and 

A 2 , , 2 



= \ [tr[C 2 ] + (1 ± 1) tr[AC 2 ]] = c 2 + ^tr[AC 2 ] , (35) 



t ^] ±tr t Z2 l ±t «]f- trM + ^). (36) 



2 V tr[Q 



d(fc 2 ) ; 

To achieve conditions CI and C2 in the 2-component BEC the effective coupling between 
the hyperfme states has to vanish, C/ 12 = 0. This can be implemented by imposing a 
particular transition rate A = —2^/p 10 p 20 Ui 2 . m addition to the fine tuning of A, the 
paramters (Uu, Pio,mi) have to be chosen so that the speed of sound simplifies to: 

c 2 _ Uu Pio _ U22 P20 _ m 2 p w U u + m 1 p 20 U 22 + UuiPwjrii + p 20 m 2 ) , . 

mi m 2 2m\m 2 

While one eigenfrequency always remains zero, Wo,i = 0, for the second phonon mode 
we get 

2 ^U 12 (p 10 m 2 + p 2 o^i)co , QC ^ 
^0,11 = j- 2 • ( 6S ) 

The mass of the modes are then defined as 

m i/u = ^0,1/11 / c o> ( 39 ) 
and thus the AM corresponds to one massless particle mi = and one massive particle, 

2 SU 12 (p w m 1 + p 20 m 2 )m 1 m 2 ^ m 2 8U 12 

[ni-zpwUn + miPwU-i-i + Uuipwmi + p 20 m 2 )} [U u + 2U 12 + U 22 \ 
propagating in the acoustic Minkowski space-time in the hydrodynamic limit. For higher 
wave numbers we obtain LIV in the form of equation (2), and the coefficients r] 2 and 774 
for the two modes are: 

ft 2 + p 20 m 2 ^0,1/11 _ / mi/n \ 2 _ ( massquasiparticies \ 2 _ . , 

772,1/11 Aci p w m 2 + p 2n m, m-imo \ M fiff / V Planck scale nf w tivfi / 



1 



-1 



7i/ii mip 10 + 7 I/n m 2 p 2 o 



(42) 

mip 20 + m 2 pi 

where 71 = 1 and 711 = m\jm 2 are dimensionless coefficients, and M e s = ^Jm\m 2 is 
defined as the analogue Planck mass. 

From the expression it is clear, that the quadratic coefficients (41) are non- universal. 
While one is always zero, 772,1 = 0, the second 772,11 depends on the interaction constant 
C/12. For C/12 ^ (C/11 + U 22 ), which is equivalent to m u <C v /m 1 m 2 , so 772,11 is suppressed. 
However, there is no further suppression — after having pulled out a factor (h~/M P {) 2 - 
for the quartic coefficients 774,1/11- These coefficients are of order one and non-universal, 
(though they can be forced to be universal, for example if 71 = 711 and the underlying 
bosons have equal masses m\ = m 2 ). 

The suppression of i] 2 , combined with the non-suppression of 774, is precisely the 
statement that the "naturalness problem" does not arise in the current model. 
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6. Summary and discussion 

We have presented an AM that can be used as a model for and motivation for 
conjectures in QGP. Low energy perturbations in a 2-component BEC with laser-induced 
transitions between the single-particle states reproduce both massive and massless quasi- 
particles in an emergent space-time. Furthermore, we have investigated higher energy- 
perturbations for a uniform condensate without background flow. This corresponds to 
particles propagating in Minkowski flat space-time with large momentum. While in the 
hydrodynamic regime the dispersion relation is LI, beyond it the dispersion relation has 
to be modified in a Lorentz violating way. 

Due to parity of the GPE we only obtained terms with even exponents of k. We 
calculated the quadratic and quartic dimensionless coefficients 772,1/11 and 774,1/11. A key 
observation is that the present model does not suffer from the naturalness problem, 
because the quadratic corrections are Planck suppressed, while at the same time the 
quartic coefficients 7/4.1/11 have no further suppression, and are actually of order unity. 
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